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Abstract 

Recent accurate measurements Q of static potentials between sour- 
ces in various representations of the gauge group SU{3) provide a 
crucial test of the QCD vacuum models and different approaches to 
confinement. The Casimir scaling of the potential observed for all 
measured distances implies strong suppression of higher cumulant con- 
tributions. The consequences for the instanton vacuum model and the 
spectrum of the QCD string are also discussed. 



1 Introduction 

The structure of QCD in its nonperturbative domain commands attention of 
theorists for many years. Confinement and chiral symmetry breaking have 
been studied both using theoretical models and lattice simulations (for review 
see i). 

However most of the models are designed to describe confinement of colour 
charge and anticharge in the fundamental representation of the gauge group 
SU (3), i.e. the area law for the simplest Wilson loop and hence linear poten- 
tial between static quark and antiquark. A supplementary and as we shall 
see below, very important information about QCD vacuum is provided by 
the investigation of interaction between static charges in higher SU(3) repre- 
sentations. Comparing static potentials for different charges one can derive 
information about field correlators in the vacuum , which is not possible to 
obtain from fundamental charges alone. 

The recent accurate measurements of the corresponding potential have 
been performed by G.Bali in M. Preliminary physical analysis of the data 
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from was presented in [Q. We present in this paper more extended in- 
vestigation and discuss new important information about the QCD vacuum 
and constraints on several QCD vacuum models. 

2 Casimir scaling of the static potential 

We define static potential between sources at the distance R in the given 
representation D as: 

Vo{R) = -Y^ra l-HW{C)), (1) 

1 — *oo ± 

where the Wilson loop W{C) for the rectangular contour C = RxT m the 
"34" plane admits the following expansion 

{W{C)) = ^Trz^Pexp {^g j A^T'^dz^ 

oo „ 

= Trz^exp / {i9Tmi)F{2)..F{n)))da{l)...da{n) (2) 

Nonabelian Stokes theorem has been used in the above expression with the 
notation F{k)da{k) = ^{xo,u^''^)E^{u^''^)T''^u'-''\ xo)da34iu^''^), where $ is 
a parallel transporter and Xq is an arbitrary point on the surface S bound 
by the contour C. The double brackets ((...)) denote irreducible Green's 
functions proportional to the unit matrix in the colour space.Q 

Since (Q) is gauge-invariant, it is convenient to make use of generalized 
contour gauge 0, which is defined by the condition ^(xo,u^''^) = 1. 

The SU{3) representations D = 3,8, 6, 15a, 10, 27, 24, 15s are character- 
ized by 3^ — 1 = 8 hermitian generators which satisfy the commutation 
relations [T",T^] = if^^T". One of the main characteristics of the represen- 
tation is an eigenvalue of quadratic Casimir operator Cd , which is defined 
according to Cj^ = T"'T"' = Cd ■ ^- Following the notations from we 
introduce the Casimir ratio do = Cd/Cf, where the fundamental Casimir 

^It makes unnecessary to write colour ordering operator in the r.h.s. of (^). Notice 
also, that the averages in refer to a single Wilson loop and do not take into account 
screening effects which add to W(C) multiloop contributions, as will be explained below 
in section 4. 
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Cp = {N^ — l)/2Nc equals to 4/3 for SU{3). The invariant trace is given by 
Tr^i = 1. 

Since a simple algebra of the rank k has exactly k primitive Casimir- 
Racah operators [0 of order mi,..,mfc, it is possible to express those of 
higher order in terms of the primitive ones. In the case of SU{3) the primitive 
Casimir operators are given by 

eg) = SabT'^T' ■ eg) = dabcT^'T'r (3) 

while the higher rank Casimir operators are defined as follows 

= d^...^)T\.r^ (4) 

where the totally symmetric tensor (i|^) .^^ on the SU{Nc) is expressed in 
terms of 6ik and dijk (see, for example, |^). 

The potential (|l|) with the definition (0) admits the following decompo- 
sition 

VniR) = dnV^'\R) + dlV^'\R) + (5) 

where the part denoted by dots contains terms, proportional to the higher 
powers of the quadratic Casimir as well as to higher Casimirs. 

The fundamental static potential contains perturbative Coulomb part, 
confining linear and constant terms 

Vd{R)=c7dR-vd-^ (6) 

The Coulomb part is now known up to two loops and is proportional to 
Cd- The "Casimir scaling hypothesis" |Ty] declares, that the confinement 
potential is also proportional to the first power of the quadratic Casimir Cd, 
i.e. all terms in the r.h.s. of (|^) are much smaller than the first one. In 
particular, for the string tensions one should get (Jd/ctf = do- 

This scaling law is in perfect agreement with the results found in 
Earlier lattice calculations of static potential between sources in higher rep- 
resentations |]lOl are in general agreement with 



To see, why this result (to be more precise - why the impressive accuracy 
of the "Casimir scaling" behaviour) is nontrivial, let us examine the colour 
structure of a few lowest averages in the expansion (^. 

The first nontrivial Gaussian cumulant in (0) is expressed through Co 
and representation-independent averages as 

Tr,,(F(l)F(2)) = {F^{1)F^{2)) = ^ (F"(1)F'^(2)), (7) 
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so Gaussian approximation satisfies " Casimir scaling law" exactly. It is worth 
being mentioned, that this fact does not depend on the actual profile of 
the potential. It could happen, that the linear potential observed in |I| 
is just some kind of intermediate distance characteristics and changes the 
profile at larger R (as it actually should happen in the quenched case for 
the representation of zero triality due to the screening of the static sources 
by dynamical gluons from the vacuum, or, in other words, due to gluelumps 
formation). The coordinate dependence of the potential is not directly related 
to the Casimir scaling, and can be analized at the distances which are small 
enough to be affected by the screening effects. [| 

Having made these general statements, let us come back to our analysis 
of the contributions to the potential from different field correlators. We turn 
to the quartic correlator and write below several possible colour structures 
for it. We introduce the following abbreviation 

^^[4]^ ^ (F"(1)T"F''(2)T^F''(3)T''F'^(4)T'^) = j'^j'brjncrjnd ^p[i]\^°-^''<^ 

where the Lorentz indices and coordinate dependence are omitted for simplic- 
ity of notation. One then gets, with some work in the last case the following 
possible structures 





~ ^ab^cd 




^p^i]^abcd 


~ ^achd 




^p[A]^abcd 


~ fadefcbe 


(FW) ~ -\N^Cn ■ i 


^p[i]^abcd 


^ facefbde 


(FW) = 


^p^i]^abcd 


~ fapmfbpnfdemfcen 





The operator Cj^ enters together with Cj^ at higher orders. Notice, that 
the terms, proportional to the square of Cr, appear in both the 66 parts 
(the first and the second strings) and higher order interaction parts (the 
last string). Mnemonically the Cd - proportional components arise from 
the diagrams where the noncompensated colour flows inside while the C|, - 
components describe the interaction of two white objects. It is seen, that the 
Casimir scaling does not mean "quasifree gluons", instead it means roughly 
speaking "quasifree white multipoles" (see discussion at the end of the paper). 

^The same is true for the general criticism of the confining potential evaluation 
on the lattice. The considered potential might be even different from linear but still 
demonstrate Casimir scaling. 
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Table 1: The Casimir-scaling and Casimir- violating string tensions and 
shifts. Based on the data from G.Bali, [hep-lat/9908021 . 



D 


ag) ■ 10^ 


Aag) ■ 10^ 


t;g) ■ 10^ 


At;g) ■ 10^ 






8 


-3.486 


1.2 


-2.513 


2.8 


0.004 


19.22 / 43 


6 


-6.428 


1.2 


0.950 


2.6 


0.007 


25.76 / 42 


15a 


-5.244 


0.55 


-0.5611 


1.1 


0.003 


39.06 / 42 


10 


-4.925 


0.50 


0.2489 


1.0 


0.003 


22.05 / 41 



Let us analyse the data from |T| quantitatively. We have already men- 
tioned, that the Coulomb potential between static sources is proportional 
to Cd up to the second loop (and possibly to all orders, this point calls for 
further study) and hence we expect contributions proportional to C|, ~ djy 
to the constant and linear terms, i.e. we rewrite (|^) as follows 

Vd{R) = doV^^\R) + ci^(t;g) + ag^i?). (8) 

and all higher contributions are omitted. Here f g"* , ag^ measure the c?!,- 
contribution of the cumulants higher than Gaussian to the constant term 
and string tension respectively. The results of the fitting of the data from |jl| 
with (j^) for some representations are shown in the Table 1. See also Fig.l, 
where the quantity {Vd^R) — doVp^R)) versus distance R is depicted. This 
figure shows the same data as the Figure 2 from the paper 

All numbers in the table 1 are dimensionless and given in lattice units. 
The author of 0] used anisotropic lattice with the spatial unit = 2AGeV . 

Several comments are in order. First of all it is seen that the Casimir 
scaling behaviour holds with very good accuracy, better than 1% in all cases 
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in the table 1 with the reasonable x^- It should be stressed, that any possible 
systematic errors which could be present in the procedure used in [|l| must 
either obey the Casimir scaling too or be very tiny, otherwise it would be 
unnatural to have the matching with such high precision. Q Nevertheless, 
the terms violating the scaling are also clearly seen. While the value of the 
constant term is found to be compatible with zero within the error bars, 
it is not the case for cr^'* . We have not found any sharp dependence of on 
the representation D, which confirms the validity of the expansion (§) and 
shows, that the omitted higher terms do not have significant effect in this 
case. Notice the negative sign of the string tension correction. In euclidean 
metric it trivially follows from the fact, that the fourth order contribution is 
proportional to {igY > while the Gaussian term is multiplied by [igY < 
for real g. 



3 Casimir scaling and instantons 

From perturbation theory it follows, as it was already mentioned that Casimir 
scaling holds up to the g^ terms. One might suspect therefore that also non- 
perturbative configurations when treated exactly, ensure the Casimir scaling. 
This is not true however for the models based on the classical solutions. As 
an example of the model which violates scaling we mention here the model 
of the dilute instanton gas. 

In the simplest SU{2) case the field strength of one instanton in the 
regular gauge is given by 



where is the position of the instanton and p is its size, chosing the Wilson 
plane to be "12" one has riai2 = ^as- The average over stochastic ensemble of 
the dilute instanton gas implies averaging over (global) color rotation of each 
instanton F12 Vt^ Fi2Vt^ averaging over instanton positions z^^'^ k = 1..N, 
where N is the total number of instantons and antiinstantons in the volume 
V and weighted averaging over instanton sizes p. The latter one is assumed to 
be performed as the last step of all calculations. To the lowest order in density 



Since we are mostly interested in the relative quantities, their actual magnituge in the 
physical units is of no prime importance for us. This is another reason why we do not 
discuss possible systematical errors of and finite volume effects. 
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(yP'^) o'^^ must consider one instanton and sum up over k,l < k < N, hence 



one can write 12 



(F(1)..F(2)) = {T'^,..Tpn ■ y J d'zF^x,, z)..F\x^, z) (10) 



The lowest order terms read 



/rp3 rp3 \ _ r . ,rp3 rp3 rp3 rp3 \ _ r 3Cj) Cq 

K-'-af^''' I3j/^ — '^a-y—y , [-L a/B-'- /Bj-'- ■yS''- Se/^ — "ae — 

Consider now the rectangular Wilson loop RxT and perform the expan- 
sion of the Wilson loop with respect to the fl- and z^ - average procedures: 

{W)n,z, = TiD (Pexptg J F^^da^^ = Tr^ exp(-A2 + A4 + ..) (11) 

The terms A2 and A4 generate the following terms in the potential: 

Vn{R) = v!^\R) + v!f\R) (12) 

where vjj^ is proportional to Co and behaves as R^ at small R, namely 
V^^\r) = 7(2)/?2/p3 where 7(2) = ^7(2) with 7(2) = ^ IQn^ For large 
distances R^ p one has Vo\r) = a^'^^R where cr^^) = -y^"^) /2p^. 
Analogously for Vo\R) one has 

(13, 

where 

'^^^''^^^ I Iw^WW^WT^+Ty ^^^^ 

Straightforward calculation gives at large distances 

T/(4) / _ ^(4) p ^(4) _ ^ ^2 1^ 3(^1) - Cp 

Vp [R) - ap R , (Jp - —p — — (15j 

while in the regime R <^ p,T ^ p one gets 

(4)_ N R' vr^ ?,Cl - Cn 

""77^ ^^^^ 
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It is clear from (115]) and (0), that linear asymptotics of vj^^ at large R 



occures rather late, for R > 7p. One concludes, that the Casimir scaling 
of the potential is violated by the term vjj\ This was interpreted in 
as the upper bound on the instanton density, for f^"* ^ 10~^G'ef it gives 
N/V ~ 0.2 fm~'^, which is much less than the instanton density typically 
used in the literature N/V ~ lfm~^. 

The original SU{3) case for the quark-ant iquark potential in the dilute 
instanton gas approximation was considered in |TB[. One has the following 
expression for the potential between static sources: 

N °r 1 /? 

V{R) = 4vr- / dpu{p)p'—- Y: (2 J + l)Fj{x) , x = f (17) 
V { d{D) f^jj 2p 

and the function Fj{x) is given by some cumbersome double integral which 
can be found in |jl3|. Here d{D) = D is the dimension of the representation D 
and sum over J = 0, |, 1, .. goes over all SU{2) multiplets for decomposition 
of the given S'[/(3) representation with the corresponding weights. One has 
E (2J + 1) = d{D) and also 

d{D) ■ Cd = ^^V^ E JiJ + 1)(2^ + 1) 

^ J&D 

At small X the functions Fj{x) ~ x^, while at large x the functions Fj{x) 
tend to J-dependent constant ||T3[| . 

Numerically one finds at small distances 

V{R) = 1.79 ■ 7i?2 . eo + 0{R^) (18) 

where 7 = ny /q°° dph'{p)p and numerical coefficients en for D = 3,8, 10 are 
given by 

63 : eg : eio = 1 : 1.87 : 3.11 

instead of Casimir scaling results 1 : 2.25 : 4.5. 

Similar situation takes place for the large distance asymptotics of the 
instanton-induced potentail. It violates Casimir scaling on the level of 20% 
(see [|^) and can be exluded by the present analysis at the level of lOa. 

So one can see the sharp contradiction between the dilute instanton gas 
model calculation for the quark-ant iquark potential and the Casimir scaling 
of this potential found on lattice. This can be understood in one of two ways. 
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Either instantons are strongly suppressed in the real(hot) QCD vacuum (as 



it was observed in 0]) while they are recovered by the cooling procedure. 
Or else instanton medium is dense and strongly differs from dilute instanton 
gas, in such a way that higher cumulant components of such collectivized in- 
stantons are suppressed. Interesting to note, that linear confinement missing 
in the dilute gas, is recovered in this case. 



4 Casimir scaling and QCD string 

There is another important consequence of the observed Casimir scaling. 
It comes from the analysis of the confinement potential as being induced 
by the QCD string. In this case one has additional contribution to the 
confining potential besides the leading linear term, which comes from the 
internal dynamics of the string, in particular, from the transverse worldsheet 
vibrations. The simplest model in this respect is the Nambu-Goto string 
which action is proportional to the area of the surface bounded by the static 
sources worldlines. It modifies the confining potential with respect to the 
classical case (nonvibrating string) as 

(tR^ (tR- — - + ... (19) 
12 R ^ ^ 

where the term — 7r/(12 ■ R) will be referred to as the String Vibration (SV) 



term Despite the Nambu-Goto string model cannot be rigorously de- 
fined in D = 4, and, in particular the expansion of the r.h.s. of (|I^) meets 
singularity at the distances i? ~ 1/ ^/a it is instructive to look whether or not 
the data support the existence of such term. It is also worth noting, that 
the dimensionless coefficient — 7r/12 is determined by the only two factors: 
target space dimension and the chosen string model. Having both factors 
fixed, it cannot be freely adjusted. Assuming an = daap, it is easy to see, 
that the Nambu-Goto induced SV term violates Casimir scaling. 

It is a nontrivial task to separate the contributions of the discussed sort 
in the confining potential as it is because these corrections are essentially 
large distance effect, where they are sub leading. But they have to become 
pronouncing in the expression (^ due to scaling violation. Namely one has 

Vn{R) - doVpiR) = [do - 1)^^ + - (20) 
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where the dots denote the terms, omitted in (0). The dashed hne on Fig.l 
corresponds to the r.h.s. of (]20|). It is seen, that Casimir scahng- violating 
SV term of the form (|I9|) is completely excluded for all measured distances.^ 
We should mention at this point, that the lattice measurements performed 
in did also demonstrate no fingerprints of the SV corrections of the form 



discussed. 

One can see different explanations of this result. First of all, nobody has 
proved up to now, that the simplest bosonic Nambu-Goto string model prop- 
erly describes the dynamics of the QCD string and theoretical background 
of ( ]I9D is not clear. Just the opposite is true - there are many resons, why 
it is not the case (see discussion in fl^)- The theory of the QCD string - 
whatever it will be - must explain the observed scaling of the potential. 



5 Screening and string breaking 

In the previous discussion the effect of string breaking in the triality zero 
representation was not taken into account. The modification of the potential 
(|^) due to this effect was considered in [|1^] in the strong coupling expan- 
sion, resulting in the expression for the adjoint Wilson loop in the large Nc 
approximation 

(WsiC)) = exp(-a8 ■ area(C)) + -i- exp(-4(T3 ■ perimeter(C)) (21) 



The appearence of the second term on the r.h.s. of ( [21| ) signals screening 
due to the lattice diagrams related to the vacuum average of the product of 
two undamental Wilson loops, while the first term is represented by the only 
one loop and has the cumulant expansion (^. We now demonstrate that a 
similar form also appears in the general background perturbation theory p2 



and estimate the corresponding gluelump masses. To this end we represent 
= Bfj^ + and use the 't Hooft identity to average separately over 5^ 
and a^, while gluon action is written as 

SiB, a) = SoiB) + S1 + S2 + S,nt (22) 



^ Notice, that even the sign of the leading Nambu-Goto string correction is opposite 
to what has actually been observed. 



10 



In this espression Si is proportional to the first power of a^, while for 5*2 one 
has 

S2 = a''(G[Br^)^y (23) 

and Sint contains higher order terms in the field a^. We integrate over Pa^, 
making perturbatve expansion in Smt, which leads to 

{W{B + a)) = J VBexp{-S{B)) [Det G[B]]-^ {W{B) + ...) (24) 

where dots stand for higher order terms in the expansion over perturbative 
fields. Note that the retained in (|2^ ) term is purely nonperturbative. Using 
the world-line Fock-Schwinger representation for DetG, one has 

Det G = exp(TrlnG) = exp I Tr J— J Vz^ exp(-K)- 

■WcAB) exp{JdT2gF)j (25) 

where the integral is to be taken with the standard Wiener measure. It is 
seen from (^5]) and (^4]) that expansion of the exponent yields expansion in 
products of Wilson loops, and the first terms are 

{WciB + a)) = {Wcm + {{WciB)WcM)) + ■■ (26) 

where the average in the last term include the average over Cz according 
to (pSj). The second term in (^) is responsible for the screening, since the 



asymptotics of it for large C is of perimeter type, rather than area law. To 
find the logarithm of (|2^) explicitly, one can use the Hamiltonian formalism, 
obtained via the transition 

Jvz exp{- J (fxL) = {x\exp{-HT)\y) (27) 

Then for the asymptotics of the second term in ( p6|) one gets 

{{WW)) = const • exp(-2MGL ■ T) (28) 

For simplified estimates we disregard the last spin term in (^) and interac- 
tion of two adjoint gluelumps with masses Mgl- To find Mgl from H one 
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can use the standard technic of einbein formahsm (see and references 
therein) to get an estimate ||18| 



Mgl ~ 1.4 GeV (29) 

This leads to an estimate of the screening distance Rq from the relation 
Vadj{Ro) = '^Mgl to be Rq ~ 1.4 Fm, which is beyond the distance where 
Casimir scaling was measured in It is interesting to note that a similar 
estimate of gluelump mass for higher D leads to the decreasing Ro{D), e.g. 
for D = 15s one gets Rq = 0.7 Fm which is not in contradiction with the 
data from |[l[]. 



6 Discussion 

The Casimir scaling behaviour of the confining potential confirmed in |jl| 
with the unprecendented precision leads to many important consequences, 
some of which have been discussed in the present paper. It is instructive 
to compare the pictures of the QCD vacuum, suggested in different models 
from the Casimir scaling point of view. 

Abehan projection language being in wide use nowadays as one of the 
most adequate for the dual Meissner scenario of confinement encounters dif- 
ficulties in explanation of the Casimir scaling. The reader is referred to the 
paper [IT3 where the question is discussed in details for the adjoint static 



charges. The observed adjoint string tension (at intermediate distances) 
arises from the interaction of diagonal abelian projected gluons with the 
part of the adjoint source doubly charged with respect to the Cartan sub- 
group. If one naively omits the corresponding Faddeev-Popov determinant 
it gives aadj = '^CFfund- It is expected, that the loop expansion of the de- 
terminant produces terms, correcting the above behaviour to the Casimir 
scaling relation. Up to the authors' knowledge, it has not yet been shown 
analytically, while there are numerical evidences from the lattice in favour 
of this possibility (see and references therein). From physical point of 
view to reproduce Casimir scaling, which is genuine nonabelian feature, one 
needs to restore the original nonabelian gauge invariance broken by hand in 
the abelian projected method. 

The now popular confining mechanism is the model of fat center vortices 



20|. While the original center vortex picture cannot explain confinement 



of the adjoint charges, the introduction of the finite thickness of the vortex 
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makes it possible, to obtain approximate Casimir ratios for the string tensions 



20| . However, with the high accuracy of the data these ratios are exluded. 

In the gauge-invariant formahsm the Casimir scahng has two im- 
portant features. First, it is the direct consequence of the Gaussian domi- 
nance hypothesis (see review p|) since gaussian correlator provides the exact 
Casimir scaling. On the other hand, it implies the cancellations of Cl, - 
proportional terms and higher ones in the cluster expansion (|^). Physically, 
it means the picture of the vacuum, made of relatively small colour dipoles 
with weak interactions between them. One can imagine two possible scenario. 
According to the first one, Casimir scaling is the consequence of Gaussian 
dominance. It this case any higher cumulant contributes to physical quan- 
tities much less than the gaussian one due to dynamical reasons. There is 
also the second possibility, when each higher term in the expansion (§) is not 
small, but their sum demonstrates strong cancellations of Casimir scaling 
violating terms. These pictures are in close correspondence to the stochastic 
versus coherent vacuum scenario 0. This set of questions certanly deserves 
further study. 

There are also several open questions of computational origin. Additional 
measurements are needed in order to clarify the validity of scaling for higher 
representations where the statistics is still rather poor. It would also be very 
interesting to establish the adjoint string breaking scale in SU (3) which could 
shed some light on the gluelumps physics. 

Needless to say, that deeper theoretical understanding of the QCD vac- 
uum structure is still required. The ability to incorporate such nontrivial 
feature as Casimir scaling is a necessary property of any reasonable confine- 
ment model. 
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Figure 1: The Casimir scaling results for the adjoint potential, based on 
the data from axes - (V8(-R) — 2.25 ■ V3(-R)); solid line - leastsquare fit 
according to (P); dashed line - scaling violating contribution from the SV 



term (|20|) . All quantities are given in the lattice units. 
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